Understanding anisotropy generated by fluctuations in heavy-ion collisions 
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Event-by-event fluctuations are central to the current understanding of ultrarelativistic heavy- 
ion collisions. In particular, fluctuations in the geometry of the early-time collision system are 
responsible for new phenomena such as triangular flow, which have solved important puzzles in 
existing data. We propose a simple model where initial fluctuations stem from independent flux tubes 
randomly distributed in the transverse plane. We calculate analytically the moments of the initial 
anisotropies (dipole asymmetry, eccentricity, triangularity), which are the sources of anisotropic flow, 
and their mutual correlations. Our analytic results are in good agreement with calculations from 
commonly-used Monte-Carlo codes, providing a simple understanding of the fluctuations contained 
in these models. Any deviation from these results in future experimental data would thus indicate the 
presence of non-trivial correlations between the initial flux tubes and/or extra sources of fluctuations 
that are not present in current models. 

PACS numbers: 25.75.Gz, 25.75.Ld, 24.10.Nz 



I. INTRODUCTION 

Fluctuations of the positions of nucleons within nu- 
clei result in event-by-event fluctuations of the density of 
matterproduced in ultrarelativistic nucleus-nucleus colli- 
sions I1h3]. These event- by-event fluctuations have been 
observed through the centrality and system-size depen- 
dence of elliptic flow Q , which is the second Fourier har- 
monic of the azimuthal distribution of outgoing particles. 
They may generate elliptic flow even in high-multiplicity 
proton-proton collisions at the LHC [5|, |6| , which might 
be responsible [3, Q for the near-side ridge observed by 
CMS 0. Furthermore, fluctuations result in new, odd 
harmonics of the azimuthal distribution of outgoing par- 
ticles, which are clearly seen experimentally through two- 
particle correlations [Ifljllll: The third harmonic, trian- 
gular flow [l2|, is largely responsible for the ridge and 
away-side structure seen in two-particle correlations and 
has recently been analyzed at RHIC [Ij] and LHC [Tjj. 
A smaller first harmonic, directed flow [15|, |la], is also 
predicted. Finally, there are non-trivial correlations be- 
tween these three harmonics |15l . 171. and they can be 
directly measured experimentally 18 1. 

Flow fluctuations are typically studied within Monte- 
Carlo models of nucleus-nucleus collisions including 
initial-state fluctuations, such as NeXSPheRIO [l9| . 
UrQMD coupled to hydrodynamics [20ll EPOS coupled 
to hydrodynamics [IJ, AMPT ^7^, and Glauber 
models coupled to hydrodynamics [24, 25\ . These com- 
plex Monte-Carlo approaches do not allow for a simple 
understanding of the key quantities driving flow fluctua- 
tions. 

In this paper, we propose a simple model of initial- 
state fluctuations and, by comparing with full Monte- 
Carlo calculations, show that it captures the essential 
physics contained in common models. We assume that 
the initial energy density profile in the transverse plane 



is the superposition of N random independent [26 
identical sources, namely: 
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where Xj are N independent random variables with a 
smooth probability distribution p{xj), and p{r) is the 
profile of a single source. The coordinate system is chosen 
such that centers of colliding nuclei are on the x axis, 
and the origin lies halfway between the two centers. For 
a collision of identical nuclei, the system then has the 
symmetries p{x, —y) = p{x, y) and p{—x, —y) = p(x,y). 
A standard choice for p{r) is a Gaussian profile [2J, |25| . 
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p{r) = poe-'-'/" 
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where a controls the width. The parameters in our model 
are iV, the number of sources, the source distribution p(x) 
and the source profile p{r). 

In Sec. ini we introduce the initial anisotropies e„ and 
their moments, which are the relevant quantities. We 
discuss to what extent e„ depends on the source profile 
p{r). In Sec. IIIIl we derive analytic expressions for the 
rms £„ and compare our results with two Monte-Carlo 
models. In Sec. IIVI our results are extended to higher 
moments of the distribution of e„ and mixed correlations 
between harmonics [l^ [13, 13 • 



II. QUANTIFYING INITIAL FLUCTUATIONS 

A. Azimuthal asymmetries of the initial 
distribution 

In a given event, the azimuthal distribution of out- 
going particles is driven by the azimuthal distribution 



of the initial density |27H29j . We define the participant 
eccentricity H £2, the triangularity £3, and the dipole 
asymmetry [l^ £1, and the corresponding orientations 
*n, by: 
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where { ■ • ■ } denotes an average over the transverse plane 
in a single event, weighted with the energy density: 



rr, s. S f{x,y)e(x,y)dxdy 
J e(x,y)dxdy 



(4) 



and (r, 0) are polar coordinates in which the pole is the 
center of the distribution, that is, {re*"^} — 0. 

The numerators in the right-hand side of E g. (l3l) are 
the cumulants introduced by Teaney and Yan [15|. The 
anisotropic flow of outgoing particles is driven by these 
initial anisotropics, in the sense that the relation w„ oc £„ 
holds approximately [lll27l-l30l| . 



£„ on the source profile is contained in the denomina- 
tor, {r"}. If the source has a finite size, {r"} increases, 
which results in a smearing of anisotropics. One expects 
that {r"} scales approximately like {r^}"/^, resulting in 
a stronger smearing for £3 than for £2, as observed in 
numerical calculations [2g. One also expects that £i/£3 
and 63/62 are largely independent of the source profile. 



C. Moments 

In practice, anisotropic flow is not analyzed in a single 
event, but in a centrality class. It is inferred from mul- 
tiparticle azimuthal correlations [3l[. To the extent that 
Vn is driven by £„, the measured azimuthal correlations 
scale like the corresponding moments of the distribution 
of £„ and their joint correlations. We introduce the no- 
tation [18|]: 



r{ni,...,nfe} = {e 
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where angular brackets denote an average over events in 
a centrality class. With this notation, the rms average of 

£„ is 



B. Dependence of e„ on source size 



,{2}^{eiy^'=s{n,-ny^^ 



(8) 



Within our simple model, one can easily study the de- 
pendence of £„ on the profile of a single source p{r). This 
is relevant to the early-time dynamics, which has a smear- 
ing effect on each source [23] and widens the source pro- 
file. 

The generating function of cumulants is W{k), 
where lla 



■,W{k} 
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Teaney and Yan's cumulants are obtained by expressing 
W{k,j;, ky) in terms of K ^ kx + iky and K = kx — iky, 
and expanding in power series of K^K''. Specifically, the 
numerators of £1, £2 and £3 in Eq. ([3]) are obtained by 
expanding W{K,K) to order KK^, K^, and K^, respec- 
tively, and {r^} is obtained by expanding W{K,K) to 
order KK. 

In our model, Eq. ([IJ, the density profile, is the con- 
volution of the profile of a single source with the distri- 
bution of sources ^ S{x — Xj). Hence 



where we have used the notation £n{2} for the rms value 

of£„ i. 

Symmetry with respect to the reaction plane (y ^■ 
—y) implies that all moments are real, therefore 
e{-ni, • • • , -Tife} = £{ni, • • • , Uk}. In addition, {x, y) -> 
{—X, —y) symmetry implies that the only nonvanishing 
moments are those such that ni + ■ ■ ■ + nk is even. Since 
the azimuthal orientation of each collision is uncontrolled 
experimentally, one can only measure rotationally sym- 
metric quantities. Thus, the only relevant moments are 
those satisfying ni + ■ ■ ■ + nk = 0- (Some of our interme- 
diate results will involve other moments, such as £{1, 1} 
or £{3,-1}.) 

Within our independent-source model, all moments 
can be calculated analytically as a systematic expansion 
in powers of 1/iV, where N is the number of sources.. 



III. RMS ASYMMETRIES 
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where p(k) is the Fourier transform of the source profile 
p(r). By symmetry, p(k) depends only on \k\'^ — KK. 
Therefore the numerators of £„ in Eq. (|31) are strictly 
insensitive to the profile p{r). The only dependence of 



The rms value of £„ is of direct relevance to analy- 
ses of anisotropic flow. Indeed, the simplest measure- 
ment of anisotropic flow is a pair correlation, (cosnAcf)), 
which is the average value over events of v'^. Assum- 
ing that Vn oc £„ on an event-by-event basis, this scales 
like (e^) = £n{2}^. In this Section, we derive analytic 
results for £^{2} within our independent-source model, 
and compare with numerical results. 



A. Analytic results 

We denote by (/(x, y)) the average value of /(a;, y) with 
the source probabihty density p{x), and we introduce the 
notation Sf = {/} — (/} for the event-by-event fluctua- 
tions. We use the complex coordinate z = x + iy. The 
asymmetry e„ is given by Eq. ([3]), where we replace re**^ 
hy z — Sz to take into account the recentering correction. 
To leading order in fluctuations, one obtains 
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(9) 



where z = x — iy. The rms value of e„ involves an average 
over events of products of i5's. Two-point averages are 
computed using the following identity, which holds for 
independent sources |2a |: 



i^f^a) = 



{f9)-{f){g) 
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We thus obtain, using the identities {z^z™) = for odd 
n — m and {z^z™) = (r"+"^ cos{{n — m)(j))) for even n — m: 
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In these equations, Ss = — (r^ cos2(/))/(r^) denotes the 
standard eccentricity. We also recall the result for £2 {2} 
which has been derived earlier 12611. 



e2{2? 
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The first term in the right-hand side is the standard ec- 
centricity, and the second term is the contribution of ec- 
centricity fluctuations to order 1/N . For a large number 
of sources, A^ ;» 1, the participant eccentricity £2 reduces 
to the standard eccentricity, while odd harmonics £1 and 
£3 vanish. 



B. Comparison with Monte-Carlo results 

We now compare analytic results derived from our 
independent-source model with results obtained using 
the mckt-vl.OO Monte-Carlo [13]. With this Monte- 
Carlo one can calculate results from both a Color- 
Glass-Condensate (CGC) inspired model — the MC- 
KLN j33| improved with running-coupling BK uninte- 
grated gluon densities [S^ , as well as a standard Monte- 
Carlo Glauber [3a]. We present only results for Pb- 
Pb collisions at 2.76 TeV per nucleon-nucleon collision. 



though results for 200 GeV Au-Au collisions agree equally 
well. In the Glauber model, each participant nucleon 
is given a weight [36[ w = 1 — x + xNcoW, where Ncoii 
is the number of binary collisions of that nucleon, and 
x = 0.18 [13. 

One input of our model is the probability distribution 
of sources in the transverse plane, p{x). For the sake 
of consistency, we assume that sources are distributed 
according to the average density profile: p{x) = {e{x)), 
where (• ■ ■ ) denotes an average over many events in a 
centrality class. We assume pointlike sources for simplic- 
ity: p(|i;|) = S{x). The last free parameter in our model 
is the number of independent sources N. One expects 
that this number scales typically like the number of par- 
ticipant nucleons in a collision. However, participants are 
not independent, but strongly correlated: for each par- 
ticipant of the projectile, there is by definition at least 
one participant from the target which is close enough in 
the transverse plane for a collision to occur. Neverthe- 
less, it is plausible that the system behaves like a set of 
N independent clusters of two or more nucleons. Again 
for simplicity, we use N = 0.45A^part for all centralities, 
though one could make the agreement with Monte-Carlo 
even better by fitting TV for each centrality class. In- 
creasing N by 20% typically decreases £3(2} and £i{2} 
by 10%, and £2(2} by less than 4%. If one uses a larger 
value of N for peripheral collisions (say, N = 0.67Vpart 
instead of A^ = 0.45Afpart), agreement with Monte-Carlo 
is significantly better for £i{2} and £3(2} but slightly 
worse for £2(2}. 
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FIG. 1. (Color online) e„{2}, with n = 1, 2, 3, versus central- 
ity. Symbols are Monte-Carlo results, lines are our analytic 
results for independent sources. For £2(2}, the dashed line is 
the full result to order 1/A'^ (Eq. (fT^ l. while the solid line is 
the standard eccentricity (first term in the right-hand side of 
Eq. dUl)). 



Fig. [T] displays a comparison between our result (|lip 



and Monte-Carlo calculations. Our analytic formulae re- 
produce the centrality dependence of e„ computed using 
Monte-Carlo KLN or Monte-Carlo GlaubeiQ- Odd har- 
monics, £i and £3, have a mild centrality dependence and 
essentially scale like l/\/iVpIrt- On the other hand, £2 
is much larger for semi-central collisions: this increase is 
driven by the almond shape of the overlap area between 
the two nuclei. For £2 , we show both the standard eccen- 
tricity, which is the leading-order term in a 1/A^ expan- 
sion, and the full expression including eccentricity fluctu- 
ations to order 1/A^ (Eq. (IT^ ). The standard eccentr icity 
is significantly larger for KLN than for Glauber [33, |40| , 
so that eccentricity fluctuations (which are roughly the 
same in both models) are a smaller relative correction. 

Our model explains why £1 < £3, a feature which 
was observed in Monte- Carlo calculations but yet unex- 
plained. This may be readily understood from Eq. fTTj) 
for central collisions, where £s = 0: only the (r^) remains 
for £3 {2}^, while £i{2}^ has an additional contribution 
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proportional to (r^)^ 



), which is negative. For sake 



of completeness, we have also carried out two other sets 
of Glauber calculations with only wounded nucleon scal- 
ing (x = 0) or binary collision scaling {x = 1). Values of 
£i/£3 turn out to be significantly smaller for wounded nu- 
cleons (between 0.3 and 0.4 for 0-50% centralities) than 
for binary collisions (between 0.6 and 0.7). This strong 
ordering is not reproduced by our analytic result, which 
gives intermediate results (between 0.5 and 0.6 in both 
cases). We do not have a simple explanation for this 
discrepancy. 



IV. HIGHER-ORDER MOMENTS 

A. Fluctuations of e„ 

The ALICE collaboration has recently measured 
V2{4:} [HI and V3{4} ^IJ]. The relative magnitude of 
z;n{4} and Vn{2} depends on event- by- event fiuctuations 
of Vn if nonfiow effects are small. Assuming that u„ is 
proportional to e„ on an event-by-event basis, fiuctua- 
tions of Vn are due to fluctuations of £„ : 



VnH} 
Vn{2} 



g»{4} 
e„{2} 
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where we have introduced the 4-cumulant £n{4} 3]. 
u„{4} is thus related to {sfj), which is e{n,n, —n, —n} 
in the notation of Eq. ([T]). Fig. [5] displays the ratio 
{£n)/i^n)'^ fo'" "■ = Ij 2, 3. The larger the ratio, the larger 
the fluctuations of £„. Note that the order of the curves 
in Fig. [2] is reversed compared to Fig. [1] smaller values 
of £,i go with larger fluctuations. For the most central 



^ We have also done comparisons with the PHOBOS Monte-Carlo 
Glauber l3al . and found that ei from this model is significantly 
smaller than predicted by our analytic formula. 
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FIG. 2. (Color online) Ratio {et)/{eiy, with n = 1, 2, 3, ver- 
sus centrality. As in Fig. [ij symbols are Monte-Carlo results. 
Full lines are leading-order analytic results. The dashed line 
for n = 3 includes next-to-leading order corrections derived in 
Appendix. Although our next-to-leading order result is valid 
for central collisions only, we compare with Monte-Carlo data 
up to 40% centrality. 



collisions, ratios are close to 2 for all n. This value cor- 
responds to Gaussian fluctuations [42|. 

For n = 2, the ratio quickly decreases with central- 
ity. This is due to the standard eccentricity Eg, which 
is zero for central collisions but dominates over eccen- 
tricity fluctuations above 10% centrality. The smaller 
eccentricity fluctuations, the closer the ratio to I. The 
decrease is stronger for the KLN model than for the 
Glauber model because of the larger Eg- Monte-Carlo 
results are compared with an analytic result using the 
expression of £2 {4} to order 1/A^ derived in [2g|. In prac- 
tice, £2 {4} — £s for all centralities. 

We now compute (£|) within our independent-source 
model. To leading order in the fluctuations, £3 is given by 
Eq. ([9]), and (£|) involves products of four S's. To leading 
order in 1/A^, average values of such products can be 
expressed in terms of two-point averages (Eq. dTU])) using 
Wick's theorem. For instance, the four-point function 
is E43| 

(SfSgShSk) = {SfSg){ShSk) + {SfSh){SgSk) + {S f S k) {S g S h) . 

(14) 
This identity gives 

£{3, 3, -3, -3} = 2£{3, -3}2 + £{3, 3}£{-3, -3}, (15) 

or, equivalently. 



(4) _,^/e{3,3} 
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£:{3, 3} can easily be calculated to order 1/A^ in the same 
way as (£3): 



£{3,3} = (£^C0S6$3) 



^, 3,2 ( (r'' COS 60} + 6£s (r ^ ) (r"* COS ' 



(17) 



Inserting Eqs. (fTTj) and ([TTl) into (fT6|) . we obtain a 
leading-order analytic result for the ratio. This analytic 
result is independent of N, and very close to 2 in prac- 
tice, as can be seen in Fig. [2l This means that the first 
term on the right-hand side is the dominant contribu- 
tion \43\. It is easy to understand why the second term 
is small: £{3, 3} involves the 6**^ Fourier harmonic of the 
initial distribution, which is of order £^. Therefore the 
last term in Eq. ([T5| is of order £^, which is very small 
in practice. 

The Monte-Carlo results in Fig. [2] show that the ratio 
(£|}/(£3)^ is slightly smaller than 2 for both Monte-Carlo 
KLN and Glauber models (a fact also implied by exper- 
imental results |44l|). which cannot be explained by our 
leading-order result Eq. ((T6)) . The next-to- leading cor- 
rection (restricted to central collisions) is derived in the 
Appendix. It is negative and scales like 1/N. As shown 
in Fig. [21 it improves agreement. 

Finally, (£f ) can be computed in the same way. The 
leading-order result is: 
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£123 = £{3, -2, -1} = (£3£2£l C0S(3$3 - 2$2 - $l)) 
£1223 = £{3, 1, -2, -2} = (£i£3£l C0S(3$3 + $1 - 4$2)) 
£13 = £{3, -1, -1, -1} = (£3£? C0S(3($3 - $l))) • (20) 

All these correlations were studied numerically in 
Ref. [l8|, with the exception of £1223 which is new. We 
derive an analytic prediction for these quantities using 
our independent-source model. To leading order in 1/iV, 
one can replace each factor of £26^^**^ in Eq. ((T]) by the 
standard eccentricity £5. The moments in Eq. (j20l) can 
then be expressed in terms of two-point moments using 
Wick's theorem: 

£{3,3,-2,-2,-2} = £M3,3} 

£{1,1, -2} =£,£{1,1} 
£{3, -2,-1} =£,£{3,-1} 

£{3,l,-2,-2} = £2£{3,l} 
£{3, -1, -1, -1} = 3£{3, -1}£{1, 1}. (21) 

To leading order in 1/A^, £{3, 3} and £{1,1} are given by 
Eqs. (|17p and ([T5)) . The remaining two-point functions 
£{3,-1} and £{3, 1} can be computed similarly: 

£{3,-l} = (£3£lCOs(3$3-$l)) 

= ]V(^(('^'c°«20)-H3£.(r2)(O 

-6£,{r2)3-3£3(r2)3-2(r2)(r4cos20) 
-|-£s (r^)(r'* cos 4(/))) 

£{3,l} = (£3£lCOs(3$3 + *l)) 
1 



(18) 
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+ 4£,(r^)(r4cos20) 
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where 



£{1,1} = (£?cos2$i) 



1 



((r^cos2</>) -8£,(r2)^ + 2£,(r2)(r^) results. We first scale £{ni, •• • ,71^} by £„,{2} • • 



iV(r3 
-4(r2)(r4cos20)-£3(r2)3). 



(19) 



This quantity is negative, which means that the dipole 
asymmetry develops mostly out of the reaction plane for 
non-central collisions. The last term in Eq. (|T8l) gives 
a positive contribution to the ratio (£i)/(£i)^, which in- 
creases up to 3 for peripheral collisions. Note that £i{4}* 
defined by Eq. ([T5)) is negative, so that £i{4} is undefined. 

The leading order result defined by Eqs. (fTS)) . ((TT|) and 
([T^ is independent of N . It is in good agreement with 
Monte-Carlo results. 

Note that within our independent-source model, the 
ratios plotted in Fig. [2] are strictly insensitive to the pro- 
file of a single source p{r). 



B. Mixed correlations 

We now study the correlations between $1, $2 and $3. 
The lowest order non-trivial correlations are 

£23 = £{3, 3, -2, -2, -2} = (£^£3 C0S(6($3 - $2))) 
£12 = £{1, 1, -2} = (£^£2 C0S(2($1 - $2))) 



These equations define our analytic results for the mo- 
ments. 

We now compare our analytic results with Monte- Carlo 

J2} 

to single out the angular correlation as in Ref. [18|. We 
compute these ratios both with the Monte-Carlo and 
with our analytic formulas. To leading order in 1/iV, we 
use the approximation £2 {2} ~ £,. One easily shows, us- 
ing Eq. (|2ip . that the resulting leading-order prediction 
for the ratios is independent of N . 

Fig. [3] displays a comparison between Monte- Carlo re- 
sults and analytic results. Our leading-order results ex- 
plain the sign, magnitude, and centrality dependence of 
all ratios. 

The scaled £23 is much smaller than unity, in agree- 
ment with earlier observations that triangularity and ec- 
centricity are uncorrelated [20| except for peripheral col- 
lisions [455. One thus expects that triangular fiow and 
elliptic fiow are essentially uncorrelated, as confirmed by 
a recent experimental analysis [1J| . Note that while data 
agree qualitatively with expectations, they do no t ag ree 
quantitatively with existing Monte-Carlo models |44| . 

We now discuss the remaining correlations, which in- 
volve $1 and are not yet measured. The scaled £13 is 
large for peripheral collisions, which means that there is 
a strong positive correlation between 3$i and 3<I'3 p7| . 
On the other hand, £12 is negative, which means that 
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FIG. 3. (Color online) Mixed correlations versus central- 
ity. From top to bottom: ei3/(£i{2}^e3{2}) (labeled 13), 
£i223/(ei{2}e2{2}2£3{2}) (labeled 1223), e23/(£2{2}'e3{2}'') 
(labeled 23), £i23/(ei{2}e2{2}£3{2}) (labeled 123) and 
ei2/(ei{2}^e2{2}) (labeled 12). Symbols are Monte-Carlo re- 
sults. Full lines are analytic results to leading-order in 1/N . 
Dashed lines for e\2 and £123 include next-to-leading order 
corrections derived in Appendix. As in Fig. [21 we display 
NLO results up to 40% centrality, although they are only 
valid for central collisions. 



$1 is more likely to be perpendicular to the participant 
plane. Because of the strong correlation between 3$i and 
3<I>3, this also exp lains why the mixed correlation £123 is 
also negative [l^ [Tj] . 

A closer look at Fig. [3] reveals that our leading-order 
results for £12 and £123 do not agree well with Monte- 
Carlo for central collisions, where both ratios are small 
but not zero, while our leading-order results go smoothly 
to zero. Agreement is improved by including the next-to- 
leading term, which is of order 1/N'^, negative, and does 
not vanish for central collisions. This term is calculated 
in the Appendix for central collisions. 



CONCLUSIONS 



sources across the transverse plane. We have considered 
pointlike sources for simplicity. The magnitudes of £3 and 
£1 show that N is roughly half the number of participant 
nucleons, which means that participant nucleons are cor- 
related pairwise. The angular correlations between $1, 
$2 and $3 are independent of N to leading order in 1/A^. 
They only depend on the distribution of sources, and are 
mostly driven by the almond shape of the overlap area 
between the two nuclei, which is responsible for the large 
elliptic flow Hi]. 

There is an ambiguity in the definition of the "initial 
state". What is computed in Monte-Carlo models (KLN 
or Glauber) is typically the distribution of energy right 
after the nuclei have passed through each other. On the 
other hand, the relevant anisotropics for collective flow 
are the anisotropics at a somewhat later time, when the 
system reaches local equilibrium. The value of £„ de- 
creases during this thermalization time [28|. We have 
shown that this decrease is solely due to the increase in 
the system size, due to the smearing of each source. £3 
depends more strongly on the system size than £2 , which 
explains why it decreases more strongly during the ther- 
malization phase. (This effect was interpreted by the 
authors of Ref. '28l| as an interference between £2 and 

Within our independent-source model, ratios such as 
(£^)/(£^)^, as well as the various angular correlations 
between event planes considered in [l8| , are independent 
of the number of sources to leading order in 1/iV (with 
the exception of {£2) / {^2y ) ■ Remarkably, these ratios are 
also unchanged through the early-time dynamics, and are 
therefore well-defined observables for initial-state fluctu- 
ations. 
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We have shown that the magnitude of initial-state 
anisotropics £„, their fluctuations and mutual correla- 
tions can be understood within a simple model where 
fluctuations stem from identical, independent sources, 
which can be viewed as "hot spots" scattered across 
the interaction region [4a, |43| . The independent-source 
model reproduces results from both CGC-inspired and 
and Glauber Monte-Carlo models. 

In our model, all information about the initial state 
is encoded in a few parameters: the number of sources 
N, the profile of a single source, and the distribution of 



Appendix A: Higher-order calculations 

In this Appendix, we present next-to-leading order 
calculations in 1/iV for a few moments. Since next- 
to-leading order calculations are considerably more in- 
volved than leading-order calculations, we assume az- 
imuthal symmetry for sake of simplicity. We compute 
three quantities for which our leading-order results van- 
ish in the limit of azimuthal symmetry, namely, £3(4}, 
£12 and £123. 

We first compute £3(4}. Taking the recentering cor- 



rection into account, the expression of £3 is 



1 



(Al) 



Treating S as the expansion parameter, we can expand 
this expression to any desired order in S. For instance, 
{{z-S^)^} = S^3-3{z^)S^-35,2S^ + 2{S^)^. We organize 
the calculation as follows. Let 



el = B + C + D 



(A2) 



where B r^ 0{S'^), C '^ 0{S^), etc. (There are no terms 
independent of S or of order S.) Hence 



and 



£3(2}^ EE {si) = (B) + (C) + (D) + 



= 2{By-{B^) 
+2{Cf - (C^) 
+2{Df - \d^) 
+i{B){C)-2{BC) 
+A{B){D)-2{BD) 
+A{C){D)-2{CD) + 



(A3) 



(A4) 



To proceed further, we make use of the expres- 
sions for {SfSg), {SfSg6h)ASfSgShSu}, (SfSgShSuSv) and 
{SfSgShSuSySw) given in [43]. (These expressions ignore 
correlations between participant positions.) Dominant 
terms in these five quantities are of 0{1/N), 0{1/N'^), 
0{l/N^), 0{1/N^) and 0{1/N^), respectively, where N 
is the number of participants. For central collisions, az- 
imuthal symmetry allows to simplify the previous expres- 
sion: 

£3{4r = 2(5)2 _ ^^2^ _ {C')+4{B){C) 

-2{BC)+A{B){D)-2{BD). (A5) 



After some algebra, we get 

"2(r6)2-(ri2) 8(r6)(r9) 



£3(4}^ 



1 

7V3 



8(r6) 



6\3 



(,3) 



3\4 



(,3) 



3\5 



(r3)e 



(A6) 

A similar result has been derived previously for £2(4}'*, 
which is also of order 1/N^ for independent sources and 
central collisions [43]. This scaling in 1/iV'^ is natural for 
the four-particle cumulant |31j . 

Other moments can be computed in a similar way. Our 
results for £123 and £12 are given in Appendix IB] 



Appendix B: Summary of results 

The expressions of £„{2}2 to order 1/N are given by 
Eqs. PT|) and P^ . The higher-order cumulant £i{4}'' 
has been derived to order l/iV^: 



£i{4r^2(£f 



2\2 



(4) 



N^{r^Y 



[(r<'cos20) - 8£,(r2)3 -I- 2£,(r2)(r4) 



-4(r2}(r4cos20}-£f(r") 



3/^2\31 



(Bl) 



The expression of £2 {4}^ was derived in [43| up to order 
1/N^ for central collisions: 



£2{4}4 ^ 2{elf {el) 

et + —^ [2ei (r^) ~ 2el {r^ cos M 



N{r^Y 



1 

iV2 



1 

7V3 



„ (r^) lr^cos26) 
8£'^-4+4£,^ -^ 



{r^f 



16£- 






16£ 



(r4)(r'* cos 2</>) 



(r2) 



,2\^ 



2{r^f (rS) 8(r4)(r6) 



(B2) 



£3{4r^2(£^)2-(£^) 
1 



N^r^) 



[{r'' 



Our result for £3{4}'* to order 1/iV'^ is 

-I- 6£,. (r2)(r'* cos 4<^) 

(B3) 



3/^2\31 



-9£^(r2) 



Ar3(r3}4 

8(r6)3 



2{r''y - {r'^) + 



(r3)2 
The 3-point mixed correlations are, to order 1/N'^, 



£12 



7V(r3; 
1 

+ 7V2 



, [(r6cos20)-8£,(r2)3 + 2£,(r2)(r4) 

-4(r2)(r4cos20)-£^(r2)3] 
4(r^)2 - (r^) 
(r2)(r3)2 
8(r2)3+4(r6)-16(r2}(r4) 



and 

£123 



iV(r3 



1 

"]V2 



(r3)2 

-[(r6cos20)+3£,(r2)(r4)-6£.(r2)3 

--3£^(r2)3 - 2 (r2)(r'' cos 20} 
-|-£s (r2)(r'^ cos 40)] 
"2(r6)-6(r2}(r4) 3{r^f - {r^) 



(B4) 



(,3) 



3\2 



(r2)(r3) 



.3\2 



(B5) 



For higher-order mixed correlations, we have only carried 
out leading-order calculations. The 4-point correlations 
considered in this paper are 



£1223 = ]yZ%2 [(''^cos40)+4£^(r2)(r4cos2<? 
-2(r2)(r'*cos40)+9£^(r2)31 



(B6) 



and 

£13 



[{r^cos2^}-8es{ry + 2es{r^){r*) 



-4{r^){r^ cos2(f,) - e^,{rY] 
X [{r^ cos2(l)) + 3es{r'^){r'^) - 6s s{rY 
^3el{ry~2{r^){r*cos2(^} 
+es{r^){r^cos4:(f))]. (B7) 



Finally, we have computed the 5-point correlation 



£23 



Nir^y 



[{r^ cos 60} + 6es (r^ ) (r^ cos 40) 



3/^2\31 



-9ei{r^) 



(B8) 



For sake of completeness, we list all other 4- and 5-point 
correlations which can be constructed out of the first 
three harmonics, and their expressions to leading order, 



obtained using Wick's theorem: 

e{2,-2,3,-3} = e2{2}2£3{2}2 

e{l,-l,2,-2} = ei{2}2£2{2}2 

£{1, -1, 3, -3} = ei{2}2e3{2}2 + e{3, l}^ + £{3, -l}^ 
£{1, 2, 3, -3, -3} = £,(2e3{2}2£{3, -1} + £{3, 1}£{3, 3}) 
£{1, 1, -1, 2, -3} = £,(2£i{2}2£{3, -1} + £{1, 1}£{3, 1}) 
£{1, 1, -2, 3, -3} = £.(2£{3, 1}£{3, -1} + £3{2}Ml, 1}) 

£{l,2,2,-2,-3} = £2{2}2£i23 
£{l,l,2,-2,-2} = £2{2f£i2 



£{l,l,l,-l,-2} = 3£i{2}^£l2. 



(B9) 



The first of these equations means that there are no cor- 
relations between the magnitudes of £2 and £3 to leading 
order. This could be tested experimentally by measuring 
the correlation between the magnitudes of V2 and V3 . 
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